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We demonstrate the possibility of three-dimensional cool- 
ing of neutral atoms by illuminating them with two counter- 
propagating laser beams of mutually orthogonal linear polar- 
ization, where one of the lasers is a speckle field, i.e. a highly 
disordered but stationary coherent light field. This configu- 
ration gives rise to atom cooling in the transverse plane via 
a Sisyphus cooling mechanism similar to the one known in 
standard two-dimensional optical lattices formed by several 
plane laser waves. However, striking differences occur in the 
spatial diffusion coefficients as well as in local properties of 
the trapped atoms. 

PACS number(s): 32.80.Pj, 42.50.Vk 



I. INTRODUCTION 

The study of speckle laser patterns, as created when 
a highly coherent light beam is transmitted through or 
reflected from an object with a surface which is rough 
on the scale of the laser wavelength, was initiated many 
years ago using electromagnetic theory JlJ and statistical 
methods Since then, this subject has raised more 

and more interest, resulting in a vast development of the 
theory (for a survey of the most recent results, see e.g. 
j|) as well as experimental achievements. These give 
rise to important applications of laser speckles in various 
fields of science, such as the speckle reduction in imagery 
||, roughness measurements in material science Q, or 
applications in bio-physics . 

The scope of the work presented here is to investigate 
the application of such speckle laser fields in the context 
of laser cooling of neutral atoms. We are especially in- 
terested in the disordered analog of the so-called optical 
lattices formed by several laser plane waves, which yield 
periodic optical potentials and have been demonstrated 
experimentally to give rise to efficient laser cooling by a 
Sisyphus- type mechanism p|-|lC|]. The basic principle of 
this cooling scheme is that the atoms loose kinetic en- 
ergy by running up potential hills from where they are 
optically pumped into lower lying potential wells. Re- 
cent experiments have also demonstrated a similar cool- 
ing scheme in the case of laser configurations forming 
quasi-periodic optical lattices fll|| , which can be viewed 
as an intermediate regime between the standard periodic 
lattices and the completely disordered patterns obtained 
from speckle light fields. 

We consider two counterpropagating laser fields of mu- 
tually orthogonal polarization, where one of the beams is 
a speckle field. According to the randomly distributed 



phase and intensity gradients of the speckle field in 
all dimensions we find three-dimensional Sisyphus cool- 
ing even with this one-dimensional beam configuration. 
Moreover, interesting transport phenomena are found in 
this case, e.g., a large difference in the spatial diffusion of 
the atoms even in parameter regimes where the steady- 
state temperatures are of the same order of magnitude 
for the longitudinal and the transverse directions. Other 
effects, such as local radiation pressure forces, arise also 
from the fact that amplitude and phase of the speckle 
field are essentially independent. 

The paper is organized as follows. In Sec. [n]we qualita- 
tively discuss our proposed setup and the basic physical 
processes which give rise to the cooling mechanism. Fur- 
thermore we discuss the dependence of the steady-state 
temperature on various system parameters. In Sec. Ill we 



present the details of the theory and outline the numer- 
ical methods, especially the semiclassical Monte-Carlo 
simulations, to obtain the results. Sec. IV is attributed 
to the analysis of the numerically obtained results such 
as temperatures, spatial diffusion coefficients and local 
effects in the light field. Finally, in Sec. fv| we discuss the 
differences of the cooling in the laser propagation direc- 
tion and in the transverse plane. 



II. QUALITATIVE DISCUSSION 

Throughout this paper we will discuss the simple situa- 
tion of a single atom with a ground state of total angular 
momentum F = 1/2 and an excited state with F 1 = 3/2 
interacting with two counterpropagating laser fields of or- 
thogonal polarization. In the case of both laser fields be- 
ing plane waves this gives rise to well-known sub-Doppler 
cooling mechanisms in one dimension (ID) |l2). Because 
of the modulation of the optical potential, trapping of 
atoms in the longitudinal direction has been predicted 
Jl3[ and observed |L4|]. In contrast, the atoms are free in 
the transverse plane. 

In the situation discussed here, one of the laser fields 
is replaced by a speckle field, i.e., by a highly disordered, 
but nevertheless stationary and coherent light field. Such 
a speckle field can be easily generated experimentally, for 
instance, by introducing a diffusor into the path of a laser 
plane wave. The resulting light field shows highly disor- 
dered intensity and phase distributions (for a discussion 
of the statistical properties of speckle fields see e.g. Q). 
An example of a computer-generated speckle field on a 
discrete spatial grid is shown in Fig. ^. 

We have investigated two different possible setups for 
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the polarizations of the counterpropagating light fields, 
namely the case of two mutually orthogonal linearly po- 
larized fields (lin_Llin configuration) and the case of two 
opposite circular polarizations (er+ — cr_ configuration). 

For the lin_Llin configuration one finds in the laser 
propagation direction essentially the usual ID Sisyphus 
cooling mechanism [ fL2| apart from the different spatial 
variation of intensity and phase of the speckle field as 
compared to a plane wave. Additionally, the properties 
of the speckle field give rise to Sisyphus cooling of the 
atom in the transverse plane, since the phase and in- 
tensity modulations in these directions lead to spatially 
varying light shifts and optical pumping rates of the Zee- 
man sublevels of the atomic ground state. Thus, the two 
counterpropagating laser beams (one plane wave and one 
speckle field) allow for three-dimensional cooling of the 
atom, in sharp contrast to the ID lin_Llin optical lattices 
used so far. 

Sisyphus cooling in the plane orthogonal to the laser 
propagation direction works analogously in the 0+ — cr_ 
configuration. However, for an f = 1/2 to F' = 3/2 
atomic transition, Sisyphus cooling in the longitudinal 
direction is expected to be much less efficient since the 
optical potentials of the ground state sublevels change on 
a much larger length scale in this direction (if the typical 
length size of the speckle in the transverse direction is 
d spi the size along the propagation (longitudinal) axis is 
d^p/^, where A is the laser wavelength). Furthermore, 
the cr+ — cr_ cooling scheme [|l2| is efficient only if F > 1 . 

In these two configurations the most striking difference 
to the configuration of counterpropagating plane waves 
will be found in the transverse plane. Hence in the re- 
mainder of this paper we will mainly restrict ourselves to 
a 2D model of cooling in the transverse directions and 
will only briefly discuss the changes induced by a model 
including the longitudinal direction in Sec. [v|. 

As an example we plot the steady-state temperature 
reached in this 2D subsystem versus the optical poten- 
tial depth HA' for fixed optical pumping rate 7' in Fig. [|, 
where A' and 7' for the speckle field are defined with re- 
spect to the mean intensity of the beam. The numerical 
method to obtain these results will be described in detail 
in Sec. 



[I] 



As our first and most important result we 
note that, as expected from the previous discussions, the 
atoms reach a steady-state temperature in the transverse 
plane. Secondly, the general behavior of this steady-state 
temperature as a function of the optical potential depth 
resembles very much the case of the standard optical lat- 
tices consisting of plane laser waves, i.e., for large val- 
ues of A' one finds a linear dependence and hence the 
steady-state kinetic energy becomes constant relative to 
the depth of the optical potentials in this limit. On the 
other hand, the temperature increases dramatically with 
small, decreasing values of A' and, finally, below a cer- 
tain threshold value of A' no steady-state temperature 
is achieved. Consequently, one finds the lowest absolute 
values for the steady-state temperature for intermediate 
values of the optical potential depth. 



According to this similarity with standard optical lat- 
tices we may roughly estimate the steady-state tempera- 
ture of the atoms in the speckle field from the formula 



k B T = 



D 



(1) 



where D is a mean momentum diffusion coefficient and 

Note that this formula 



a mean friction coefficient |15 



only holds in the absence of atomic localization, which 
is not the case here as we will see later, but nevertheless 
gives a useful order of magnitude for the temperature. 
The friction coefficient will be of the order of 



A 



A 



(2) 



where A is the detuning of the laser from the atomic reso- 
nance frequency, T is the natural linewidth of the excited 
state, and d sp is the mean distance of neighboring inten- 
sity maxima along a ID cut through the speckle field, 
i.e. the typical length scale of the speckle field. Eq. (Q) 
is obtained from the well-known expression for standard 
lattices [ ^2p^ | where we have only replaced the factors of 
A by d sp , which accounts for the different typical length 
scale. 

The momentum diffusion coefficient can be guessed as 



D = D 



dip 



D s 



Irlr [ -^ A' 
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-,+in 2 k 2 , (3) 



which contains the diffusion Ddi P due to the fluctuating 
dipole force and the diffusion D se due to the recoil of the 
spontaneously emitted photons. Ddip is approximately 
given by the square of a typical force, Uk{\/ d sp )(A /T) , 
times a typical time, I/7', and D se has the same form as 
for a standard lattice. Hence the estimated temperature 
(FtT) reads 



k B T 
HA' 



T 



(4) 



This formula exactly predicts the qualitative behavior of 
the temperature in Fig. |^, i.e. the linear increase for large 
values of A' and a rapid increase for very small values. It 
also shows that the lowest value of T (= T min ) is achieved 
for AY7' = d sp /X and that /csT mi „ = 2h r y'd sp /X. For a 
fixed optical pumping rate 7', the minimum temperature 
is thus expected to increase linearly with the speckle size. 

Finally, one notes from Fig. || that, although both the 
lin_Llin and the <r + — <r_ configuration give rise to cool- 
ing, the lin_Llin setup is more efficient and hence gives 
rise to lower temperatures. For the chosen parameters 
the difference in temperature between these two setups 
is about a factor of two in the vicinity of the threshold 
value for the optical potential depth and is smaller for 
larger values of A'. 

Another important difference between the laser cool- 
ing inside a speckle field as discussed here and standard 
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optical lattices created by several plane waves lies in the 
tunability of the typical length scale. For our setup the 
mean speckle grain size d sp can be changed continuously 
by changing the position of the diffusor which creates the 
speckle field out of a plane laser wave. 

As an example, Fig. || shows the dependence of the 
steady-state temperature on d sp for the lin_Llin as well 
as for the <r + — <r_ configuration. One clearly sees 
an increase of the temperature for larger speckle sizes. 
Eqs. (||)-(||) explain these feature since we note that the 
friction force decreases as l/d^ p , but only the dipole 
diffusion term Ddi P has the same dependence, whereas 
the sponaneous emission term D se remains constant. 
Hence the latter contribution of the total diffusion starts 
to dominate for large speckle sizes (more precisely for 
d sp /\ 3> A 1 /Y), and hence the temperature is expected 
to increase quadratically with the speckle size. However, 
for the parameters chosen in Fig. || (A'/t' = 15 ) the tem- 
perature increases only by a factor of two, if the speckle 
grain size is increased from about A to 12A. But as we 
will discuss later in Sec. IV, the increase, for instance, of 
the spatial diffusion coefficient, which is equal to D/a 2 
in a certain range of parameters fl7|| , is much larger. 

Again the qualitative behavior of the er + — er_ config- 
uration is the same as that of the lin_Llin configuration. 
Thus in the remainder of the paper we will only focus on 
the latter situation, since it is more appropriate for a 3D 
cooling mechanism as discussed previously. 



III. THEORETICAL MODEL 

In this section we give a more detailed discussion of the 
mathematical model and the numerical methods used to 
obtain the results presented in this paper. 

A. Speckle fields 

The first step for all the numerical treatments is to 
create a speckle field. This can be easily implemented 
on a computer following the procedure described in the 
work of Huntley , which we will briefly outline in the 
following. 

Let us denote the electric field at the diffusor by 
Ei(x, y) and the electric field at the plane chosen for the 
numerical simulations by E 2 (x, y). The two fields are re- 
lated according to the Huygens-Fresnel principle in the 
Fresnel approximation by 



where I denotes the distance between the diffusor and the 
observation plane. Hence, if we define 



E 2 {x,y) = Yi exp 



x / dx 1 dy 1 E 1 (xi 1 yi)exp "^(^l+yi) 



x exp 



'.2tt, 

1 AT Xl +VVl > 



(5) 



Ei(x,y) = E 1 {x,y)exp \-i—(a? + y») 
E' 2 (x, y) = E 2 (x, y) exp \i ^(x 2 + y 2 ) 



(0) 
(7) 



these two quantities are related by a Fourier transform. 

Following H the real and imaginary part of E' 2 (x 7 y) at 
any point (x, y) are independent Gaussian random vari- 
ables. Thus the numerical construction of a speckle field 
starts by filling the real and imaginary parts of E' 2 (x,y) 
on a discrete spatial grid of N x N points with random 
numbers from a Gaussian distribution of zero mean and 
unit standard deviation. The grid size N must be chosen 
in such a way that it contains a reasonable large number 
of speckle grains, but nevertheless the discretized field 
must be smooth enough in order not to introduce large 
numerical errors. In practice we have used values of N 
between 64 and 256 for our calculations. 

The effect of the finite size of the diffusor is imple- 
mented by transforming E 2 (x, y) into E'^x, y), multiply- 
ing the latter quantity by a window function W(x,y) 
which assumes unity inside the diffusor and vanishes 
otherwise, and transforming the result back into a final 

E , 2 M = F- 1 (WF(E' 2 )), (8) 
where T denotes the two-dimensional Fourier transform. 



B. Semiclassical model of atomic dynamics 

As in most of the previous theoretical works on 
the Sisyphus cooling mechanism, including semiclassical 
fll2"| , [l6" 19 1 as well as quantum treatments j2Tjj2l[], we con- 
sider the simple case of an atom with a ground state 
of angular momentum F — 1/2 and an excited state of 
angular momentum F' = 3/2. Furthermore we restrict 
ourselves to a 2D model and to the case of low atomic 
saturation, where we can adiabatically eliminate the ex- 
cited state of the atom. The time evolution of the atomic 
density operator restricted to the ground state manifold 
is then governed by the master equation 



P = --AH,p] + Lp, 
with the Hamiltonian 

H = p 2 /2m + hA'V(x) 
and the decay and recycling term 

Lp = i/2[-V{x)p-pV{x) 

+ 2W d 2 qN{q)e-^Bl{x)pB CJ { 



(9) 



(10) 



iqx 



(11) 
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where the sum goes over the polarization and the integral 
over the wave vector of the spontaneously emitted photon 
projected into 2D. The optical potential depth is KA', 
the optical pumping rate is 7', and the atomic transition 
operators B a (x) and the optical potential operator V(x) 
are defined as 



4, (12) 

A a = J2( F '^'\^^ F ^)\F',m')(F,m\, (13) 

m' ,m 

V{x) = Y,B a {x)Bl{i), (14) 



where E a (x) gives the spatial dependence of the a- 
polarized laser light and where we have made use of the 
Clebsch-Gordan coefficients in Eq. (|l3|). 

In order to derive a semiclassical theory suitable for 
Monte-Carlo simulations we rewrite the master equation 
(0) in the Wigner representation defined by 



W(x,p,t) = 



1 



(2k) 



d 2 u(x + u/2\p(t)\x - u/2)e- ipu . 

(15) 



Note that W(x,p, t) is still an operator in the Hilbert 
space of the internal atomic degrees of freedom. However, 
for our specific choice of an atomic F = 1/2 to F' = 3/2 
transition and the laser polarizations always lying within 
the same plane, no coherences between the ground-state 
sublevels build up. Hence the Wigner operator W(x,p, t) 
remains diagonal, 

W(x,p,t) = W+(x,p,t)\m = 1/2) (m = 1/2 1 

+ W-(x,p,t)\m = -l/2)(m= -1/2|. (16) 

For these diagonal terms we obtain the Fokker-Planck 
equations 



W± + — diW± 
m 



-J±t w ±+Jt± w t 
+ F± ± d Pz W ± + F^ ± d Pi W T 
+ D% ± d Pi d Pj W ± + D% ± d Pi d Pj Wzp, (17) 

where i, j — x,y and the sum over i, j must be performed. 
The full expressions for the jump rates 7± T between the 
ground-state sublevels, the force coefficients F±± and the 
diffusion coefficients D±_j_ are given in the appendix. 
Finally, we use semiclassical Monte-Carlo simulations 
to obtain numerical solutions of Eqs. p7|), where one 
follows the trajectories of many particles with internal 
states |+) and |— ). The jump rates between these two 

states are given by 7^ and 7 |_, respectively. Between 

two jumps the particles evolve according to the force F+ + 
(resp. F^__) acting on them and in addition receive ran- 
dom kicks which are chosen in such a way as to simulate 
the effect of the diffusion terms -D+ + (resp. D % 1_). 



Averaging over a set of particles and over time yields 
all the required expectation values such as temperatures, 
mean local velocities, position distributions or spatial dif- 
fusion coefficients. We will discuss the most important 



results obtained in that way in Sec. IV 



C. Estimate of final temperature 

Before turning to the numerical results obtained from 
solving the equations of motion presented in the previous 
subsection by Monte-Carlo simulations, we will analyti- 
cally derive a rough estimate for the steady-state tem- 
perature T obtained as the ratio of the mean momentum 
diffusion coefficient D over the mean friction coefficient 



k B T = — , 

a 



(18) 



where the bar over D and a denotes averaging over the 
internal atomic state and over position. For simplicity 
we will restrict the following calculations only to the x 
direction. 

In order to calculate the friction coefficient a we must 
find the stationary solution W(x, p) of Eq. ( p~T|) up to first 
order in the atomic velocity v = p/m, i.e., we expand 
W(x,p) by 



W(x,p) = W°(x) + ^-W^x) + 



(19) 



and insert this into Eq. (|l7]) . Since we are considering an 
atom moving with constant velocity we may neglect the 
force and diffusion terms on the right-hand side of the 
equation and thus obtain the results 



wl = 1 - w°_ = — ^± . 

7-+ + 7+- 
1 , d x Wl 



1-+ + 1+- 



(20) 
(21) 



The friction coefficient, i.e. the term of first order in the 
velocity v of the force, averaged over the internal atomic 
state is then given by 



a = Fl.Wi+F^_WL 



(22) 



Finally this must be averaged over the position within 
the speckle field. Since F±± as well as W± can be ex- 
pressed in terms of the speckle electric field amplitude 
and its spatial derivatives (see appendix), we need to 
know the expectation values of products of these quanti- 
ties when averaged over position. Considering the com- 
plete randomness of the speckle field we assume that all 
averages over such products vanish except for 



E+El = E-E*_ = 1, 



(23) 



assuming that the speckle field and the counterpropagat- 
ing plane wave have the same average intensity, and 
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d x E+d x E* + =d x E_d x E*_ 
2 V d m 



(24) 



where again d sp is the mean speckle grain size. The 
right-hand side of Eq. ( 24 ) is obtained under the assump- 
tion that the required quantity for the speckle field is 
the same as for a periodic electric field with the same 
typical length scale, i.e. for an electric field given by 
E(x) = cos[kx(X/d sp )]. Applying these assumptions to 
Eq. ( p2| ) yields the friction coefficient 



a = nk —— — — 

4r \d sp 



(25) 



The averaged diffusion coefficient has two different 
components. One of these is obtained by averaging the 
diffusion coefficients D±°± over the internal atomic states 
and over position similarly as done above for the friction 
coefficient. The second contribution to the total diffu- 
sion coefficient arises from the change of the dipole force, 
if the atom changes its internal state. This dipole dif- 
fusion term can be derived easily following the lines of 
Ref. jl5). We will omit all the calculational steps here 
and only give the final result for the total momentum 
diffusion coefficient 

D I [A\ 2 / A \ 2 11 / A \ 2 5 ,„„. 



h 2 k 2 -y' 



36 \d 



36' 



where the first term is the dipole term and the second 
and third terms are due to the momentum diffusion of 
an atom in a definit, internal state, i.e., coming from 

u ±±- 

Thus we obtain for the estimated steady-state temper- 
ature (lisl) 



k B T 2 , n /r\ 2 5 /r\ 2 ( d. 



hA' 



3 27 V A 



+ 



27 V A 



(27) 




We see that the expressions derived in Eqs. ( |25| ) 
essentially the same as the intuitive ones of Sec. 

It should be emphasized again that all quantities de- 
rived in this subsection should be considered as crude es- 
timates since they are performed in ID and rely on some 
very rough approximations. However, the qualitative be- 
havior of the exact solutions is predicted correctly and 
thus these expressions provide a lot of physical insight. 



IV. NUMERICAL RESULTS 

A. Steady-state temperatures 

First, let us continue the discussion of the numeri- 
cally found steady-state temperature depending on vari- 
ous system parameters. 



In Fig. [| we plot the steady-state temperature ver- 
sus the detuning A for fixed optical potential depth A', 
i.e., for varying optical pumping rate 7'. The two lower 
curves correspond to the same optical potential depth 
but to different speckle grain size d sp . In agreement with 
Eq. ( p7|) the temperature increases for smaller detunings 
A, i.e. larger pumping rates 7'. For very large detunings, 
corresponding to small values of 7', the temperature as- 
sumes a constant value, which is the same for different 
speckle sizes as long as the optical potential depth A' is 
the same. On the other hand for fixed speckle size but 
different values of A' the temperatures achieved in the 
limit of small values of 7' differ (cf. also Fig. ^|). 

Fig. U shows the steady-state temperature versus the 
optical potential depth A' for fixed optical pumping rate 
7', similar to Fig. ^ but for different parameters. Again 
we find a drastic increase of the temperature towards 
small values of A'. In this parameter limit the total mo- 
mentum diffusion ( p6| ) reduces to its contribution stem- 
ming from the random recoil of the spontaneously emit- 
ted photons and hence becomes constant, but simulta- 
neously the Sisyphus cooling vanishes, i.e., the friction 
force (|2^) approaches zero. Hence the diffusion starts 
to predominate the cooling effect and the temperature 
increases. 

In the opposite limit of large values of A' the total 
momentum diffusion is dominated by the dipole diffusion 
term of Eq. (26) and hence the temperature increases 
linearly with A'. Note on Fig. || that the steady-state 
temperature finally becomes independent of the optical 
pumping rate 7' and of the speckle size d sp in this limit, 



in agreement with Eq. (E 



B. Local properties 

In this subsection we will discuss some of the localiza- 
tion properties of the cold atoms. To this end we give 
a contour plot of the steady-state atomic density p in 
Fig. ^|(a), where the calculations have been performed 
for the speckle field depicted in Fig. [I]. Comparing these 
two figures one can see easily that the atomic density 
is strongly correlated with the speckle field intensity I. 
This statistical dependence can be demonstrated more 
clearly by calculating the covariance, defined as 



cov(p,7) = <p/)/vVK7 2 ). 



(28) 



We find that this quantity assumes values close to its 
maximum value of one, e.g. for Fig. ^(a) cov(p, I) — 0.92. 
This means that the intensity distribution of the speckle 
field is efficiently mapped onto the atomic density distri- 
bution. 

However, if one looks closer at Figs, [lj ||(a), one notes 
some discrepancies, for instance some of the relatively 
shallow optical potential wells at the top of the figures are 
surprisingly strongly populated. Although these features 
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cannot be explained quantitatively, we will describe some 
of the physical mechanisms in the following. 

Basically these local differences are caused by the fact 
that the field amplitude and the phase of the speckle 
field are essentially uncorrelated. Whereas in usual op- 
tical lattices the light amplitude and phase are always 
correlated, which e.g. in a 3D lin_Llin setup yields that 
the optical potential minima coincide with places of pure 
circular polarization ||, this is not the case in our setup. 
Consequently the local polarization of the total light field 
at places of maximum local speckle field intensity can be 
linear as well as circular, but in general will be an ar- 
bitrary elliptical polarization. Hence only speckle grains 
(regions of high speckle field intensity) of mainly circular 
polarization contribute to the Sisyphus-type cooling, but 
speckle grains of dominating linear polarization do not. 

It should be emphasized that in speckle grains of 
mainly circular polarization local cooling can be found, 
that is, even within a single well an atom can be cooled 
by optical pumping processes between the widely differ- 
ing optical potentials of the two ground-state sublevels. 
This is in contrast to the plane-waves ID lin_Llin config- 
uration, where an atom must travel across several poten- 
tial wells in order to get cooled (nonlocal cooling). Thus, 
the atom is more likely to be cooled and trapped within 
circularly polarized speckle grains. 

Another reason for the different steady-state atomic 
densities in potential wells of the same depths is formed 
by the radiation pressure force. As already mentioned 
above, the light field amplitude and phase of the speckle 
field are mutually independent and thus also the phase 
gradient can differ strongly at the bottom of the poten- 
tial wells. Hence, within some potential wells the atom 
experiences a relatively strong radiation pressure force 
and is thus pushed away. As an example, Fig. ^|(b) 
shows again a contour plot of the speckle field intensity 
(same as Fig. [j]), together with the vector field of the 
local radiation pressure force. This shows, e.g., a sig- 
nificant radiation pressure in the x direction in the high 
intensity speckle grains at the center of the figure. Cor- 
respondingly, we have numerically found nonvanishing 
mean atomic velocities in these regions, which means that 
there exists a stationary flow of atoms following roughly 
the direction of the local radiation pressure force. This 
might suggest that a description of speckled lattices in 
terms of fluid mechanics formalism may be interesting to 
investigate, but, however, this is out of the scope of this 
paper. 

Finally we will briefly discuss the atomic density dis- 
tribution P(p), i.e., the probability of finding a certain 
density. A histogram plot for the speckle field intensity 
P(I) and the atomic density is given in Fig. 0. 

The numerically obtained speckle field intensity dis- 
tribution, Fig. 0(a), follows an exponential law, as theo- 
retically discussed in Q, with equal mean and standard 
deviation, i.e., 



As can be seen from Fig. f?](b) the atomic density does not 
follow the same law, since there is always a background 
of unbound atoms and hence the probability of finding 
a density below a certain threshold vanishes. But above 
this threshold the density distribution function closely 
resembles the speckle intensity distribution, which again 
indicates a strong correlation between the light field in- 
tensity and the atomic density. 



C. Spatial diffusion 

Another important quantity for characterizing the 
properties of the cooled and trapped atoms is the spatial 
diffusion coefficient, which roughly quantifies the trans- 
fer of atoms between several potential wells and hence of 
the spreading of an initially small atomic cloud. 

Experimentally, as well as in our Monte-Carlo simula- 
tions, the spatial diffusion coefficient D s is obtained from 
the spreading of an initial atomic cloud whose variance 
follows a linear law in the long-time limit, 



[Ar(t)] 2 = [Ar(0)] 2 + 2D s t. 



(30) 



An estimate for the value of D s can be derived by using 



D, 



D 

7^2 



(31) 



for the spatial diffusion coefficient for a Brownian mo- 
tion which in certain limits is a good approximation for 
Sisyphus coolin g [fl7| . Applying our findings for the fric- 
tion coefficient ([25|) and the momentum diffusion ( p6| ) we 
obtain 
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P(//<J))=exp(-//(/». 



(29) 



We can now compare this with the spatial diffusion 
coefficients obtained by semiclassical Monte-Carlo simu- 
lations as depicted in Fi g, as a function of the detuning 
A. It follows from Eq. @ that for ^'/A' = T/A < I 
the first term will dominate and hence the spatial diffu- 
sion coefficient scales linearly with 7', i.e. like 1/A if A' 
is kept constant. This is in qualitative agreement with 
the curves of Fig. || for not too large values of the detun- 
ing A. However, Eq. ([32]) only holds if the atom moves 
only a fraction of the speckle size d between two optical 
pumping processes juj, which explains the deviation of 
the numerically obtained curves from the 1/A law for 
large values of A. 

Note that for A ^> T and within its range of validity 
Eq. (|3^) predicts an increase of the spatial diffusion co- 
efficient with the square of the speckle grain size. This 
has been nicely confirmed by our numerical simulations. 



G 



For an actual experiment, where the speckle size is typi- 
cally of the order of tens of wavelengths, this means that 
the spatial diffusion of the atoms in the speckle field can 
easily be orders of magnitude larger than for a standard 
optical lattice, even if the steady-state temperatures ( p7| ) 
are comparable. 

V. INFLUENCE OF THE LONGITUDINAL 
DIRECTION 

Compared to the transverse directions, which have 
been discussed so far, the physical properties of the sys- 
tem along the laser propagation direction are completely 
different. First, it should be emphasized again that the 
size of the speckle grains in this direction is generally 
much larger, i.e., by a factor on the order of d sp /X. But 
superimposed on this relatively large length scale there is 
a standard ID lin_Llin configuration with a typical length 
scale of half an optical wavelength, such that the system 
must be characterized in this direction by two widely dif- 
ferent length scales. Second, this superimposed lin_Llin 
configuration leads to efficient cooling in the longitudi- 
nal direction at regions in space where the speckle field 
intensity approximately equals the intensity of the coun- 
terpropagating plane wave. However, at locations where 
the two beam intensities differ it also leads to a non- 
balanced radiation pressure force driving the atoms away 
from these locations. Note that the situation is similar 
to what can be found in quasi-periodic lattices [fll| . 

Hence, the full 3D situation is much more intricate 
than the transverse 2D scheme. Also numerically a com- 
plete 3D treatment would be very demanding. Instead we 
performed some 2D simulations of the cooling dynamics 
including the longitudinal and one transverse direction, 
which should give some relevant information about the 
actual behavior of the system in 3D. 

In Fig. ^ we depict the temperatures in the longitu- 
dinal direction and in the transverse direction obtained 
from such 2D Monte-Carlo simulations. We see that the 
transverse temperatures are always higher than the lon- 
gitudinal ones. Even if there exists a strong radiation 
pressure in the longitudinal direction in regions of differ- 
ing speckle and plane wave intensity, this effect is over- 
compensated by the much more efficient Sisyphus cooling 
in this direction. 

Secondly, let us remark that the transverse temper- 
atures are larger for this type of Monte-Carlo simula- 
tions than those obtained from simulations including two 
transverse directions. Hence the inclusion of the longitu- 
dinal direction tends to reduce the efficiency of the trans- 
verse cooling. This can be understood by the following 
arguments. Since the transverse cooling relies on a Sisy- 
phus effect, it is more efficient for deeper potential wells 
because in this case a single quantum jump between the 
two adiabatic potentials reduces the atomic energy by a 
larger amount. Thus the cooling mechanism is most effec- 
tive in regions of maximum speckle field intensity. But 



in the longitudinal direction, these regions correspond 
exactly to the regions of maximum radiation pressure, 
and hence the atom is efficiently pushed away. For the 
transverse cooling this means that the atom avoids the 
deepest potential wells, which consequently reduces the 
cooling efficiency and therefore increases the steady-state 
temperature. 

For the spatial diffusion coefficient we found a much 
larger difference between the longitudinal and the trans- 
verse direction than for the temperature. This agrees 
well with our previous results that the leading term in 
the expression of the temperature ( p7| ) is independent of 
the speckle size whereas the spatial diffusion (^) scales 
with the square of the speckle size. Hence the difference 
in the typical length scale between the longitudinal direc- 
tion (A) and the transverse direction (d sp ) yields widely 
different spatial diffusion. Again we note that the val- 
ues for the transverse spatial diffusion obtained by the 
Monte-Carlo simulations with the longitudinal and one 
transverse direction lies above the value obtained from 
simulations restricted to the transverse plane, according 
to the reduced cooling efficiency as discussed above. 

VI. CONCLUSIONS 

In this work we have demonstrated that 3D laser cool- 
ing of neutral atoms can be achieved by a ID setup, if 
one of the two counterpropagating laser beams of a plane- 
waves ID optical molasses |14j is replaced by a speckle 
laser field. The cooling mechanism in this case relies on 
a Sisyphus effect similar to the one known in usual (pe- 
riodic) lattices p"2| in one, two, and three dimensions. 
We have calculated steady-state temperatures similar to 
those obtained for periodic lattices for large red detun- 
ings of the lasers from the atomic resonance. In this far- 
detuned limit the steady-state temperatures also become 
independent of the mean speckle grain size. 

In contrast to usual lattices we have found various local 
effects, such as local cooling and nonvanishing radiation 
pressure forces which are essentially decoupled from the 
local optical potential depth. Furthermore, according to 
the speckle field the spatial diffusion of the atoms is in- 
creased by a factor of approximately the square of the 
typical speckle grain size. 

Thus for an experimental realization of this scheme 
two points should be emphasized. First, the experimen- 
tal setup is relatively simple since a ID configuration is 
already sufficient to yield 3D sub-Dopplcr cooling. On 
the other hand, the obtained speckle grain size is typi- 
cally much larger than an optical wavelength which gives 
rise to longer cooling times and a huge increase of the spa- 
tial diffusion. Thus, the lifetime of the trapped atomic 
cloud will be much shorter than in usual lattices and the 
achievable atomic densities much lower. However, first 
experimental results |22| have already confirmed most of 
the theoretical predictions presented in this work and rea- 
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sonable agreement with the numerically obtained values 
for, e.g., the steady-state temperatures has been found. 

Although such investigations and possible extensions of 
the scheme, e.g., to the study of speckled dark lattices, 
already present interesting results for laser cooling, we 
think that a main interest lies in the possibility of creat- 
ing a disordered sample of cold atoms with well controlled 
statistical characteristics. Because of the growing inter- 
est in the study of disordered materials, the presented 
system may thus be useful in this domain. 
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APPENDIX A: THE FOKKER-PLANCK 
EQUATION 

In the following we give the exact expressions for the 
coefficients which appear in the Fokker-Planck equation 
( |l7j ) as obtained from the Wigner transform of the mas- 
ter equation (^) (the derivation is straightforward but 
lengthy and will thus be omitted). The jump rates read 



+ c.c. 



7 ± T=g7'l^| 2 , 

the force coefficients are 
Fi ± = A'|(^ ± )4 + i(^)4 

+7'~{(^£±)i4 + £(w4 

F^ ± =i l -{{d l E ± )El-c.c), 
where i = x,y, and the diffusion coefficients are 
DT± = l'\{^d x E ± ){d x El) + l(d x d x E T )El 



(Al) 



-~E^{d x d x E%) + -(d x E T )(d x E^)} 



1 



(A4) 



n yy 
u ±± 



7' 



'\{mjE±)(d y Ei) + 2 -{d v d v E T )E^ 
^E T (d y d y E^ + ^(d y E T )(d y E^)} 

-7'^|£ ± 4^| 2 , 
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n xy 
JJ ±± 



n yx 
u ±± 



,1 



+ ^{d x E T ){d y E%)+c.c), 



D l 4± = -i^^d^Ei+E^mEi) 

-2(^)^4)} + 7 '^|£ ± | 2 , 
-r'^{(d x d v E ± )Et 
(d x E ± )(d v El)+c.c.y 



n xy 



n yx 



(A7) 



(A8) 



For the derivation of the diffusion coefficients we have 
simplified the spontaneous emission pattern by assuming 
that fluorescence photons are only emitted along the x, 
y, and z axes Q . 

Two other points are worth a comment here concern- 
ing the application of these diffusion constants and forces 
in t he sem iclassical Monte-Carlo simulations outlined in 
Sec. Ill B| . First, it turns out that for all realistic choices 



of the system parameters the cross terms and D l ^ ± 
between the two ground-state sublevels can be neglected 
as compared with the terms F±± and D±±. Second, 
at some positions in space the diffusion coefficients as- 
sume negative values. This property indicates a purely 
quantum feature of the system, i.e., it shows that there 
exist positions in space where the atomic wavefunction 
collapses rather than spreads out p3|| . However, this fea- 
ture cannotbe mimicked by the semiclassical simulations, 
where negative diffusion makes no sense. Hence, when- 
ever this situation occurs in our numerical simulations 
we set the diffusion coefficients equal to zero to avoid 
this problem. 



-{2(d x E ± ){d v El) + -{d x d y E T )E% [8] 
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FIG. 4. Temperature (2D) vs detuning A for fixed poten- 
tial depth. Solid curve: d sp = 1.9 A, A' = 200^, dashed 
curve: d sp = 5.7 A, A' = 200 lur, dotted curve: d sp — 1.9 A, 
A' = 1000ojr. 

FIG. 5. Temperature (2D) vs optical potential depth for 
fixed optical pumping rate. Solid curve: d sp — 3.8 A, 
7' = Gujr, dashed curve: d sp — 7.6 A, 7' = 6 lur, dotted 
curve: d sp — 3.8 A, 7' = 20 lur. 

FIG. 6. (a) Contour plot of the atomic density for 
A' = 200 lur and 7' = 13.33 lur. (b) Contour plot of the 
speckle field intensity and vector field of the local radiation 
pressure force. The size of the spatial region is (10A) 2 . 

FIG. 7. (a) Histogram plot of a numerically obtained 
speckle field intensity. Theoretical analyses predict an ex- 
ponential behavior (smooth curve), (b) Histogram plot of the 
atomic density for A' = 200 lur, 7' = 13.33 lur, d sp — 1.9 A. 



FIG. 8. Spatial diffusion vs detuning A for fixed poten- 
tial depth. Solid curve: d sp = 1.9 A, A' = 200 lur; dotted 
curve: d sp — 1.9 A, A' = 1000 lur; dashed curve: d sp = 5.7 A, 
A' = 200 lur. (For the corresponding steady-state tempera- 
tures see Fig. 0). 



FIG. 9. Steady-state temperatures in the transverse direc- 
tion (solid line) and in the longitudinal direction (dashed) for 
fixed optical pumping rate 7' = 20 lur versus optical potential 
depth. The speckle grain size is d 3p — 1.8 A transversally and 
d's P = 7.2 A longitudinally. 



FIG. 1. Contour plot of the intensity of a numerically cre- 
ated speckle field on a grid of 64 x 64 points. 



FIG. 2. Temperature T vs optical potential depth HA' for 
the lin_Llin configuration (solid line) and the cr+ — <r_ con- 
figuration (dashed) for fixed optical pumping rate 7' = 6 lur. 
The average speckle g rain size is d sp — 3.8 A. 



FIG. 3. Temperature vs speckle grain size for the lin_Llin 
configuration (solid line) and for the a+ — <j_ configuration 
(dashed) for A' = 200 ljr and 7' = 13.33 lur. 
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